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Abstract. We study the performance of a multidimensional matched filter as a 
follow-up module of the coherent method recently developed by two of us for the 
detection of gravitational wave bursts by spherical resonant detectors. We have tested 
this strategy on the same set of injections used for the coherent method and found 
that the matched filter sensibly improves the determination of relevant parameters as 
the arrival time, amplitude, central frequency and arrival direction of the signal. The 
matched filter also improves the false alarm rate, reducing it roughly by a factor of 3. 
The hierarchical structure of the whole analysis pipeline allows to obtain these results 
without a significant increase of the computation time. 
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1. Introduction 

Whether the first direct gravitational wave (GW) detection wiU be made during the 
forthcoming joint LIGO- Virgo [1] run, or rather wiU be claimed by the next generation 
of interferometers, the intriguing era of gravitational wave astronomy is in any case 
expected to begin. 

In this perspective, spherical resonant detectors, having the capability of locating on 
their own the gravitational wave source in the sky, may represent an important resource 
for this arising discipline. Contrarily to an interferometer, a sphere is a multichannel 
detector capable in principle to determine all the four parameters characterizing a GW 
(two angles for the arrival direction, and the amplitudes of the two polarizations); 
in order to deal with this distinguishing feature, several data analysis issues must be 
however tackled (see j2]- [12] for an incomplete but elucidating historical perspective). 

In [13] and jH] a coherent burst search algorithm has been proposed, that fully 
exploits the multichannel capabilities of the sphere by transforming the detector's 
outputs into a five-dimensional time series corresponding to the spheroidal components 
of the gravitational wave. This time series is in turn processed by a coherent, wavelet- 
based module [15j (henceforth CWS, as Coherent Waveburst for the Sphere) generating 
a list of candidate burst events endowed with all the relevant parameters (arrival time, 
arrival direction, amplitude, etc..) and is also able, to some extent, to reject non- 
transverse excitations (which are surely not due to GW's), thus reducing the false alarm 
rate by a factor of about 10. 

The reason of using wavelets rather than the traditional matched filter(MF)-ing 
technique is two-fold. MF-based searches are computationally expensive, as they require 
the exploration of a parameter space which is made extremely large by the presence of 
several channels (for example, the multichannel MF is direction dependent, diflFerently 
from the traditional single channel one). Moreover wavelet-based searches, contrarily 
to a MF one, do not require any detailed prior knowledge of a signal, a very welcome 
feature for an unmodeled burst search such as the one developed in [131 EJ • 
On the other hand, as the use of the MF is known to give the optimal signal-to-noise ratio 
(SNR) p^, the wavelets fiexibility has a cost in term of accuracy in the determination 
of the burst parameters. 

To this respect, the present work can be considered as a natural evolution of [T3l[T4] 
in that it implements a multidimensional coherent MF as a follow-up of the method cited 
above. This hierarchical strategy is crucial in reducing the parameters space because 
for each trigger only a little region is explored around the values indicated by the CWS 
module, thus making the whole procedure reasonably fast. 

The plan of the paper is the following: in section [2] the whole pipeline is briefiy 
described, showing in particular how the MF module is interfaced to the previous part 
of the analysis chain; then in section [3] the improved method is tested against the 
same mock data sets (software injections of signals on a background endowed with both 
Gaussian and non-Gaussian noise, generated on the basis of a numerical simulation of 
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the relevant components of the detector), thus showing the improvements obtained in 
the parameters estimation. FinaUy, section [4] contains our conclusions. 

2. Description of the analysis pipeline 

2.1. From the transducers^ outputs to the event candidates 

The first part of the pipeline, which is exhaustively discussed in and fT3], will be just 
briefiy summarized here, leaving aside many technical details. The first step consists 
in converting the six transducers' outputs Ik into five channels corresponding to the 5 
quadrupolar modes of the GW tensor hij: 

Ik{t)^hN{t), (1) 

with 

= E yini.M , Yi'l = J2 y!nL^'^' ' (2) 

A^=0,4 ij=x,y,z 

(2) 

where ttin = 0, Ic, Is, 2c, 2s labels the components of the real spherical harmonics l^iv 
(defined and normalized as in [3j), y^l^^ij denotes the Cartesian component of the tensor 
spherical harmonics and is the versor of the arrival direction of the GW. 

This conversion allows to extract the information about the GW from the detector's 
outputs (or, equivalently, to express its noisy content in term of its "GW-equivalent" ) 
and is possible once a mathematical model for the detector is available. In the case at 
hand, the model consists of a set of coupled oscillators describing the dynamics of the 
relevant sphere vibrational modes, of the transducers, and of the electrical circuit that 
are at the core of the readout devices (see [14j for details). 

The next step is to combine the Fourier transforms of the five Hn channels into a 
single channel H: 

hM{uj)^H{uj). (3) 

This process is performed by extracting from the dataset the spectral density matrix, 
defined from 

S%^\uj)5{uo - J) =< hN{uj)h^^;{u) > , (4) 
and then by building the following quantity: 

Hicu) = {h-h)s=J2 ^!v(^) • ■ hN'iou) . (5) 

This quantity is then fed into the CWS algorithm, which is a modified version of the 
Coherent Waveburst algorithm developed for LIGO ([I7]-[20j), that has been adapted 
to our computational needs (a standard MacPro machine with two 3GHz processors) 
and to our detector configuration. This algorithm provides a list of triggers containing 
the various event features (central time of events, amplitude, arrival direction, central 
spectral frequency and spectral frequency spread). 
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The trigger list is finally further refined through an independent determination of 
the arrival direction (the so-called determinant method), and a trigger is vetoed when 
the two reconstructed directions are not compatible. To summarize: 

H{ljj) — > {triggers} — > {event candidates} . (6) 

2.2. Follow-up of the events via multidimensional matched filtering 

If the time dependence, arrival direction and polarization of the GW signal were known, 
that is if the expression for the five functions h^^^{uj) were available, the optimal detection 
strategy would be to perform a multidimensional matched filter. As explained in [6j, 
the multichannel MF consists in analyzing the following scalar channel: 

HMF{uj) = {h''' ■h)s, (7) 
and is the strategy which gives the optimal SNR 

S^R^ f(h-'-h)sA^ (8) 



for the given signal /i^^^. The integrals are done on both positive and negative 
frequencies. However, in a real experimental situation /i^^^ is not known a priori and 
to employ efficiently a MF one should explore a huge number of points in a parameter 
space describing the temporal profile of the signal, as well as the arrival direction and 
polarization. Thus the use of a MF in such a context is computationally very demanding. 

The CWS method exposed in the previous paragraph has been shown to find 
candidate events with a good efficiency /reliability ratio, without need to know the form 
of the signal. This allows to keep under control the computation time but of course it 
comes at a price for what concerns parameter estimation. 

The improvement we are introducing in this work consists in performing a 
multidimensional MF only in the time segments where candidates events have been 
selected by the CWS, and to restrict the analysis to a small portion of parameter space 
centered around the parameters provided by the CWS triggers. 

Specifically, we consider a two parameter family of sine-Gaussian shapes to describe 
the time profile of the signal 

hf{t\to, A, /o, e, 4>) = /ioe-(*-*°)'/(2^') sin [27r/o(t - to)] 4>) , (9) 

impinging from the generic direction characterized by polar angles ^, 0. Since we have 
verified that the unmodeled search alone already gives a good estimation of the frequency 
spread l/(27rA) (an estimation which is not improved by the MF procedure), we kept it 
as fixed, in such a way that, for each candidate event we are basically exploring a small 
portion of the 4-dimensional parameter space {to, /o, ^5 0}- 

Such a space is explored through a maximum-seeking algorithm, which searches for 
the point characterized by the maximum SNR. 

We maximized first on the arrival time, as the peak of the SNR is very sharp 
against this parameter, then on the central frequency and finally over the two angles of 
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/i„,(Hz-i/2) 


#GW 


#non-GW 


10-20 


48 


8 


7 ■ 10-21 


33 


5 


5 ■ 10-21 


30 


7 


3 • 10-21 


24 


9 



Table 1. Detection efficiency for injected GW signals and non-GW disturbances. 
Adapted from [14j . 



the arrival direction together, keeping at each stage the other parameters fixed at the 
value indicated by the CWS triggers. This is a simplified choice of course, with respect 
to building a full parameter template bank, for instance, but it allows to obtain good 
results while keeping a low computational burden. 

Once found, the maximum-SNR point in the parameter space defines the improved 
parameters of the candidate event, and if they happen to be too much different from 
the ones given by the unmodeled search, the event is discarded as spurious. 

This means that the MF follow-up procedure, as well as improving the events' 
parameters, acts as a further veto against non-GW triggers, i.e. against signals which 
are not compatible with the transverse-traceless nature of GW's. Such false triggers can 
come from Gaussian noise or from non-Gaussian disturbancies. The CWS method alone 
is already optimized to reduce the latter's, and we will show that the application of the 
MF follow-up improves the performance with respect to the rejection of false triggers 
due to Gaussian noise. 

3. Results 

We have tested the MF method on the triggers obtained in [H], i.e. 50 GW-injections, 
(25 for each polarization) for each of the four hrss amplitudes: 3, 5, 7, 10 x 10~^^Hz~^/^ 
and also on an equal number of non-GW injections (signals whose power has been 
distributed randomly among the five sphere channels) with the same amplitude. The 
injected signals are sine-Gaussian of the same type of the matching function given in 
Q, this allows us to assess the precision in parameter estimation irrespectively of a 
mismatch between the filter and the signal in the data. An eventual mismatch could of 
course worsen the estimations, but we do not investigate here on that. 

For the sake of clarity, we have reported in table [T] the detection efficiency for both 
GW and non GW-injections obtained in [Hj; these numbers, which show that the CWS 
alone already reduces by a factor of 5-10 the number of non-Gaussian false triggers, have 
not changed in the present analysis. 
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3.1. Efficiency vs. false alarm rate 

The GW detection efficiency alone is not a satisfactory measure of the goodness of the 
analysis method, as it must be completed by informations about the false alarm rate. 
A Receiver Operation Characteristic (ROC) curve has been constructed in order to 
find the set of analysis parameters (the most important being in our case the maximal 
allowed distance between the two direction reconstruction methods, Sthr) that provides 
the optimal efficiency vs. false alarm performance. 

In figure [T] we report the ROC curves both for the CWS unmodeled search and 
after the application of the MF, for different injection strengths. The false alarm rate 
is obtained from the analysis of simulated data containing Gaussian noise only. The 
figure shows that, at a given level of the efficiency, the MF method allows to reduce 
considerably the false alarm rate (by roughly a factor of 3 for 6thr = 0.2 radians, which 
the same working point chosen in [Mj). 




} . . . — . — ^ 0,0 

0,001 0,010 

False alarm rate [Hz] 



Figure 1. ROC curves (detection efficiency vs. false alarm rate) for the GW- 
injections, obtained by varying the maximum allowed distance between the direction 
reconstruction methods. Broken lines show the ROC obtained in with the 
unmatched method, varying the distance between the CWS and the determinant 
method; soHd Hues show the results after the application of the matched filter 
and are obtained varying the distance between CWS and MF methods. Lines of 
increasing thickness correspond to injections sets of increasing amplitude, hrss = 
{3,5,7,10} X lO-^^Hz-^/^ 
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3.2. Parameters estimation 



3.2.1. Arrival time The MF improves by an order of magnitude the arrival time estima- 
tion, as can be seen in figure [2| The results are summarized in table [2} where we plotted 
the average and the standard deviation of the quantities tews — Unj and tMF — t 
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Figure 2. Histograms showing the arrival time reconstruction precision for the 
unmodeled (light) and matched-filtered (black) search, at different amphtudes. 



/i,,,(Hz-i/2) 


Atu (msec) 


AtMF (msec) 




a 




a 


10-20 


43 


28 


-1.2 


2.1 


7 • 10-21 


43 


25 


-0.48 


2.9 


5 • 10-21 


52 


28 


-0.24 


3.4 


3 • 10-21 


48 


29 


-0.40 


4.7 



Table 2. Average /i and standard deviation a of the difference (in mihiseconds) 
between the arrival time determination and the actual injected one, for the unmodeled 
search (U) and for the matched fiher. The fact that < a in the MF case means 
that our estimator is unbiased, and in such case the precision of the method is given 
by a. 
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3.2.2. Central frequency Also the estimation of /o benefits of an one order of 
magnitude improvement with the MF method, see table [3| More specifically, 
most of the improvement comes from the correction of a few very imprecise 
estimations given by the CWS, as can be seen in the histograms of figure [3 . 
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Figure 3. The four histograms show the distribution of the quantities /q^^ — /q^^ 
(light) and /q^^ — /q^"^ at different injected amplitudes. 



/i,,,(Hz-V2) 


A/o^ (Hz) 


A/tf^ (Hz) 




a 
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10-20 


0.28 


4.6 


0.08 


0.28 


7 • 10-21 


0.60 


4.7 


0.08 


0.32 


5 • 10-21 


0.70 


4.1 


0.10 


0.27 


3 • 10-21 


0.97 


3.6 


0.06 


0.54 



Table 3. Average /i and standard deviation a of A/o (in Hz) for the unmodeled 
search (U) and for the matched filter. The fact that < <j means that our estimator 
is unbiased, and in such case the precision of the method is given by a. 
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3.2.3. Arrival direction In fig. [4] the direction reconstruction of the unmodeled 
and of the MF search are plotted. In this case the improvement, although clearly 
visible, is not as spectacular as in the previous cases, and appears to be more 
marked for (j) than for 9; this last feature might depend on the chosen injection 
direction, and further study is necessary to understand the behavior for generic arrival 
directions. Quantitative results for the angular distance Ss are summarized in tab. [4} 
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Figure 4. Arrival direction reconstruction by the unmodeled search (crosses) and 
by the MF (triangles) for the four given signal amplitudes. The circle represents the 
injection direction. 
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5 • 10-21 
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3 • 10-21 
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Table 4. Average angular distance (in radians) between the direction determination 
and the actual injection direction, for the unmodeled search (U) and for the matched 
filter. Differently from the previous cases, the angular distance is clearly a biased 
parameter (because it is not distributed around zero, by definition), so the accuracy 
of the estimation is well described by the average. 
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n 


a 


10-20 


2.9 ■ 10-3 


0.057 


7 • 10-21 


3.5 ■ 10-3 


0.066 


5 • 10-21 


2.0 ■ 10-3 


0.061 


3 • 10-21 


2.5 ■ 10-3 


0.086 



Table 5. Average /i and standard deviation cr of {h^^^ — h^Q^)/h^Q^ . The fact that 
< <j means that our estimator is unbiased, and in such case the precision of the 
method is given by a. 

3.2.4' Amplitude The CWS method have not been exphcitly cahbrated to measure 
the amphtude of the signal. However, it is stiU interesting to see that with the use of 
the Matched Filtering the relative precision in the determination is a few percent s, as 
reported in in table [5} 

4. Conclusions 

Matched filtering techniques are widely used in gravitational wave analysis. However 
they can have a high computational cost, because of the possibly large volume of 
parameter space required by a comprehensive search. The case of spherical detector 
is even worse as its multimode feature makes it highly sensitive to the arrival direction, 
adding the source position in the sky as further parameter. 

Here we have shown how, in the case of burst searches, combining an unmodeled search 
with a matched one can exploit the goodies of both while taming the baddies of both. 
The unmodeled makes no assumption on the type of signal, but gives a non-optimal 
estimation of the parameters of the burst events triggered. On the other hand the 
matched filtered search is optimal once the shape of the signal is known. Here we 
use the output triggers of the unmodeled search to massively restrict the search in 
the parameter of the best matching template; more specifically we used a 4-dimensional 
parameter space (arrival time, central frequency and source angles in the sky) to present 
an example of matched filtering search. Our search is very much simplified with respect 
to a general one as the parameter space is not fully explored. 

To make a comprehensive analysis one would need template banks and systematic 
scanning of the parameter space, but even with this primitive attempt we have shown 
how the first estimate made by the unmodeled search can be improved by a very light 
matched filter search. 
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